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I.  Introduciion 


TtKre  are  roughly  four  clas&es  of  lechnique^  (o  predici  aerodynamic  tbrcea  and  moments  on  a 
pro^ciile  in  atmospheric  flight:  empirical  methods,  wind  tunnel  testing,  computational  fluid 
dynamics  simulation,  and  spark  range  testing.  Empincal  techniques  aenodynamically  describe 
the  projectile  with  a  set  of  geometric  properties  (diameter,  number  of  fins,  nose  type,  nose  radius, 
etc  )  and  catalog  aerodynamic  coefficients  of  many  different  projectiles  as  a  function  of  these 
features.  The  database  of  aerodynamic  coefficients  as  a  function  of  projectile  teanires  is 
Typically  obtained  from  wind  tunnel  or  spark  range  tests.  This  data  is  fit  to  multivariable 
equations  to  create  genenc  models  for  aerodynamic  coefficients  as  a  function  of  these  basic 
pro^ctile  geometric  properties.  Examples  of  this  approach  to  projectile  aerodynamic  coefficient 
estimation  include  Missile  DATCOM,  PRODAS^  and  AP9k  {1-6).  The  advantage  of  this 
technique  is  that  it  is  a  general  method  applicable  to  any  projectile.  However,  it  is  the  least 
accurate  method  of  the  four  methods  just  mentioned,  particularly  for  new  configurations  that  fall 
outside  the  realm  of  projectiles  used  to  form  the  basic  aerodynamic  database.  The  empincal 
method  has  been  found  very  useful  in  conceptual  design  of  projectiles  where  rapid  and 
inexpensive  estimates  of  aerodynamic  coefficients  are  needed. 

In  wind  tunnel  testing,  a  specific  projectile  is  mounted  in  a  wind  tunnel  at  vanous  angles  of 
attack  with  aerodynamic  forces  and  moments  measured  at  vanous  Mach  numbers  using  a  sting 
balance.  Wind  tunnel  testing  has  the  obvious  advantage  of  being  based  on  direct  measurement  of 
aerodynamic  forces  and  moments  on  the  projectile.  It  is  also  relatively  easy  to  change  the  wind 
Tunnel  model  to  allow  detailed  parametnc  effects  to  be  investigated.  The  main  disadvantage  to 
wind  tunnel  testing  is  that  it  requires  a  wind  tunnel,  which  is  modestly  expensive.  Furthermore, 
dynamic  derivatives,  such  as  pitch  and  roll  damping,  as  well  as,  Magnus  force  and  moment 
coefficients,  are  difficult  lo  obtain  in  a  wind  tunnel  and  require  a  complex  physical  wind  tunnel 
model  Wind  tunnel  testing  is  often  used  during  projectile  development  programs  to  converge  on 
fine  details  of  the  aerodynamic  design  of  the  shell  (7,5). 

In  computational  fluid  dynamics  (CFD)  simulation,  the  fundamental  fluid  dynamic  equations  are 
numerically  solved  for  a  specific  configuration  The  most  sophisticated  computer  codes  are 
capable  of  unsteady  time  accurate  computations  using  the  Navier-Stokes  equations.  Examples  of 
These  tools  include,  for  example,  CFD*-*,  Fluent,  and  Overflow-D.  Over  the  past  couple  of 
decades,  tremendous  strides  have  been  made  in  the  application  of  CFD  for  prediction  of 
aerodynamic  loads  on  air  vehicles,  including  projectiles.  This  numerical  technique  is  based  on 
first  principles  and  does  not  involve  physical  testing.  It  is  a  general  method  that  is  valid  for  any 
projectile  configuration.  However,  CFD  is  computationally  expensive  and  requires  powerful 
computers  to  obtain  results  m  a  reasonably  timely  manner  (9-2?).  In  spark  range  aerodynamic 
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lestmg,  a  [irojeciile  is  fired  through  an  enciosed  bmidmg  Ai  a  di&creie  number  of^tomts  during 
the  flight  of  (he  [irojeciile  (<30)  the  state  of  the  projecitle  measured  using  sfark  shadowgraphs 
(23-27)  The  projectile  state  data  is  subsequently  fit  to  a  rigid  6  degree-of- freedom  projectile 
model  using  the  aerodynamic  coefTicients  as  (he  ruling  parameters  (2f(-30)  Spark  range 
aerodynamic  testing  is  considered  the  gold  standard  for  projectile  aerodynamic  coefllcient 
estimation  It  is  the  most  accurate  method  for  obtaining  aerodynamic  data  on  a  specific 
pro^ciile  configuration  It  is  usually  the  most  expensive  aliemaiive,  requires  a  spark  range 
facility,  and  is  only  valid  for  (he  specific  projectile  configuration  tested. 

Various  researchers  have  used  CFD  to  estimate  aerodynamic  coefficient  eaiimation  of 
projectiles  Early  work  focused  on  Euler  solvers  applied  to  steady  flow  problems,  while  more 
recent  work  has  solved  the  Reynolds  averaged  Navier-Stokes  equations  and  large  eddy 
simulation  Navier-Stokes  equations  for  both  steady  and  unsteady  conditions  [9-22).  For 
example,  to  predict  pitch  damping,  Weinacht  prescribed  projectile  motion  to  mimic  a  typical 
pitch  damping  wind  tunnel  test  in  a  CFD  simulation  to  estimate  the  difTerent  components  of  the 
pitch  damping  coefficient  of  a  fin-stabilized  prc^ectile  (Si)  Excellent  agreement  between 
computed  and  measured  pitch  damping  was  attained.  Algorithm  and  computing  advances  have 
also  led  to  coupling  of  CFD  codes  to  projectile  ngid  body  dynamics  codea  for  simulation  of  free 
(light  motion  of  a  projectile  in  a  time-accurate  manner.  Aerodynamic  forces  and  moments  are 
computed  with  (he  computational  fluid  dynamics  solver,  while  the  free  flight  motion  of  the 
proj  ectile  is  computed  by  integrating  the  rigid  body  dynamic  ( RBDl  equations  of  motion.  The 
ability  to  simulate  the  flight  of  a  projectile  using  first  principles  has  led  to  (he  notion  of  ^'virtual 
fly  outs”  where  the  simulation  tools  are  used  to  replicate  a  spark  range  test.  Along  these  lines, 
Sahu  achieved  excellent  agreement  between  spark  range  measurements  and  a  coupled  CFD/RBD 
approach  for  a  tinned  stabiliaed  projectile  (32)  Projectile  position  and  orientation  at  down-range 
locations  consistent  with  a  spark  range  test  were  extracted  from  the  output  of  the  CFD/RBD 
software  to  compute  aerodynamic  coeflicients.  Standard  range  reduction  software  was  utilized 
for  this  purpose  where  good  agreement  was  obtained  in  comparison  to  example  spark  range 
results. 

While  coupled  CFD  'RBD  simulation  is  now  capable  of  replicating  time  accurate  projectile 
motion,  computing  time  for  this  type  of  analysis  is  exceedingly  high  and  does  not  currently 
represent  a  practical  method  for  typical  flight  dynamic  analysis,  such  as  impact  point  statistics 
(circular  error  probable)  computation  where  thousands  of  fly-outs  are  required  Furthermore, 
this  type  of  analysis  does  not  allow  the  same  level  ofunderstanding  of  the  inherent  underlying 
dynamics  of  the  system  that  ngid  body  dynamic  analysis  using  aerodynamic  coefficients  yields. 
However,  the  coupled  CFD/RBD  approach  does  offer  an  idea)  way  to  rapidly  compute  (he 
aerodynamic  coefTicients  needed  for  ngid  6  degree-of-freedom  simulation.  Dunng  a  time- 
accurate  CFD/RBD  simulation,  aerodynamic  forces  and  moments,  and  (he  full  rigid  body  state 
vector  of  the  projectile  are  generated  at  each  time  step  in  the  simulation.  This  means  that 
aerodynamic  forces,  aerodynamic  moments,  position  of  the  mass  center,  body  orientation, 
translational  velocity,  and  angular  velocity  of  the  projectile  are  all  known  at  the  same  time 


2 


insuni  With  time-synchronized  air  load  and  stale  vector  information,  the  aerodynamic 
coefTicienis  can  be  estimated  with  a  simple  fitting  procedure.  This  report  creates  a  method  to 
evidently  generate  a  complete  aerodynamic  model  for  a  pro;ectile  m  atmospheric  flight  using  n 
short  time  histones  at  m  different  Mach  numbers  with  an  industry  standard  time-accurate 
CFD/l^BD  simulation  The  technique  is  exercised  on  example  CFD^  RBD  data  for  a  small  fin- 
stabilized  projectile.  The  technique  is  further  exercised  for  a  fot-  and  spin-stabilized  projectile 
using  simulated  data  from  a  standard  trajectory  code.  Parametric  trade  studies  investigating  the 
number  of  time  snippets  and  the  length  of  each  time  snippet  to  obtain  accurate  aerodynamic 
coefficients  are  reported. 


2 .  Projec  (ile  CFD/KBD  Simu  I  a  lion 


The  projectile  CFDyRBD  algorithm  employed  here  combines  a  rigid  6  degree-of-freedom 
projectile  flight  dynamic  model  with  a  three-dimensional  (3-D),  time-accurate  CFD  simulation. 
The  RBD  dynamic  equations  are  integrated  forward  in  ome,  where  aerodynamic  forces  and 
moments  that  drive  motion  of  the  projectile  are  computed  using  the  CFD  algorithm 

The  RBD  projeciile  model  allows  for  3  translation  degrees  of  heedom  and  3  rotation  degrees  of 
h^dom.  As  shown  in  figures  I  and  2,  the  /  frame  is  attached  to  the  ground  while  the  d  frame 
IS  fixed  to  the  projectile,  with  the  fg  axis  pointing  out  the  nose  of  the  projectile  and  the  Jg  and  Kg 
unit  vectors  forming  a  right -handed  triad.  The  projecdle  state  vector  is  comprised  of  the  menial 
position  components  of  the  projectile  mass  center  2 ),  the  standard  aerospace  sequence  Euler 

angles  ( ^,  (ir ),  the  body  fl’ame  components  of  the  projectile  mass  center  velocity  ( u,  v,  w ),  and 
the  body  frame  components  of  the  projectile  angular  velocity  vector  ( p,q,r )  Both  the 
translauonal  and  rotational  dynamic  equations  are  expressed  m  the  projeciile  body  reference  frame 
The  standard  rigid  projectile  and  body  frame  equations  of  motion  are  given  by  equations  1-4. 
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Figure  I  RefcrctKe  frami;  and  poaiiion  definitions 
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Note  ihai  ihe  total  applied  force  components  [X,Y,Z)mA  moment  component  iL,M  yN) 
contain  contributions  from  weight  and  aerodynamics.  The  aerodynamic  portion  of  the  applied 
loads  in  ei^uations  3  and  4  is  computed  using  the  CFD  simulation  and  passed  to  the  rigid  body 
dynamic  simulation 


On  the  other  hand,  the  CFD  flow  equations  are  integrated  forward  in  time,  where  the  motion  of 
the  projectile  that  drives  flow  dynamics  is  computed  using  the  RBD  algorithm.  The  complete  set 
of  3-D  lime-dependent  Navier-Siokes  equations  are  solved  in  a  time-accurate  manner  for 
simulation  of  free  flight.  The  commercially  available  code,  CFD*  *17-20,  is  used  for  the  dme- 
accurate,  unsteady  CFD  simulations.  The  basic  numerical  framework  in  the  code  contains 
unified-grid,  unified-phystcs,  and  unihed-computing  features.  The  3-D,  time -dependent 
Reynolds-averaged  Navier-Stokes  (RANS)  equations  are  solved  using  the  following  finite 
volume  equation. 
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(5> 


where  If  i&  ihe  vector  of  conservative  vanabies,  Fand  (!>  are  the  inviscid  and  viscous  flux 
vectors,  respeciively,  H  is  die  vector  of  source  lerms,  y  is  ihe  cell  volume,  and  ^  is  the  surface 
area  of  the  cell  face.  A  second-order  discretizaiion  is  used  for  the  flow  variables  and  ihe 
turbulent  viscosity  equation.  The  turbulence  closure  Is  based  on  topology-parameter-^e 
formulations.  Two-equaiion,  higher-order  RANS  turbulence  models  are  used  for  the 
computation  of  turbulent  flows  These  models  are  ideally  suited  to  unsunictured  book  keeping 
and  massively  parallel  processing  due  lo  dieir  independence  from  constraints  related  lo  the 
placement  of*  boundaries  and/or  zonal  interfaces. 
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A  dual  iime-stepping  approach  i&  to  sncegraie  (he  flow  equations  lo  achieve  ihe  desired  tune 
accuracy.  The  first  is  an  *'oiJter"  or  global  (and  physical)  time  step  that  corresponds  lo  (he  time 
discreitaaiion  of  the  physical  lime  vanauon  term.  This  lime  siep  can  be  chosen  directly  by  ihe 
user  and  is  typically  set  to  a  value  lo  represent  1/100  of  the  period  ot*  oscillation  expected  or 
forced  in  (he  transient  flow.  It  is  also  applied  to  every  cell  and  is  not  spatially  varying  An 
artificial  or  "inner"  or  "local"  time  vanauon  term  is  added  to  the  basic  physical  equations  This 
time  step  and  corresponding  "inner-iteration"  strategy  is  chosen  to  help  satisfy  (he  physical 
transient  equations  to  (he  desired  degree  For  (he  inner  iterations,  (he  time  step  is  allowed  to  vary 
spatially.  Also,  relaxation  with  muliigrtd  (algebraic)  acceleration  is  employed  to  reduce  the 
residues  of  the  physical  transient  equations,  (i  is  found  that  an  order  of  magnitude  reduction  in 
the  residues  is  usually  sufiicteni  to  produce  a  good  transient  iteration. 

The  projectile  in  the  coupled  CFD/l^BD  simulation,  along  with  its  grid,  moves  and  rotates  as  the 
pro^ctile  flies  downrange.  Grid  velocity  is  assigned  to  each  mesh  point.  For  a  spinning  and 
yawing  projectile,  the  grid  speeds  are  assigned  as  if  the  gnd  is  attached  to  the  projectile  and 
spinning  and  yawing  with  it 

In  order  to  properly  initialize  the  CFD  simulation,  two  modes  of  operation  for  (heCFD  code  are 
utilized,  namely,  an  uncoupled  and  a  coupled  mode.  The  uncoupled  mode  is  used  to  initialize 
the  CFD  flow  solution,  while  the  coupled  mode  represents  the  time-accurate  coupled  CFD;TIBD 
solution.  In  the  uncoupled  mode,  the  rigid  body  dynamics  are  specified.  The  uncoupled  mode 
begins  with  a  computation  performed  in  "steady-state  mode,”  with  the  grid  velocities  prescribed 
to  account  for  the  proper  initial  position  (  Xq  .  >'•>  Zd  )•  orientation  ( 4$) ,  dig ,  ),  and  transbiional 

velocity  (Ug,  Vg,  H*g )  components  of  (he  complete  set  of  initial  conditions  to  be  prescribed.  After 
(he  quasi-steady-siate  solution  is  converged,  the  initial  spin  rate  ( 1  is  included  and  a  new 
quasi-steady-staie  solution  is  obtained  A  sufficient  number  of  time  steps  are  performed  so  that 
(he  angular  oneniaiion  for  (he  spin  axis  corresponds  to  the  prescribed  initial  conditions.  This 
steady-state  flow  solution  is  the  starting  point  for  the  coupled  solution.  For  (he  coupled  solution, 
(he  mesh  is  translated  back  to  (he  desired  initial  position  )  and  (he  remaining  angular 

veloary  initial  conditions  (^•,rn )  ^  added.  In  (he  coupled  mode,  the  aerodynamic  forces 
and  moments  are  passed  to  the  RBD  simuboon  which  propagates  (he  rigid  state  of  the  projectile 
forward  in  time. 


3.  Flighl  Dynamics  Projectile  Aerodynamic  Model 


The  applied  loads  in  equations  3  and  4  are  expressed  in  the  body  reference  frame  and  split  Into 
contributions  due  to  weight  and  body  aerodynamic  force  as  shown  in  the  following; 
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TtK  air  loads  can  be  further  split  into  a  steady  air  loads  component  that  acts  at  the  center  of 
pressure  and  a  Magnus  air  loads  componeni  that  acts  at  the  center  of  Magnus.  The  terms 
containing  consiiniie  the  Magnus  air  loads  component,  while  the  terms  containing 

define  the  loads  acting  at  the  center  of  pressure.  The  externally  applied  moment 
about  the  projectile  mass  center  is  composed  of  an  unsteady  aerodynamic  moment,  along  with 
terms,  due  to  the  foci  that  the  center  of  pressure  and  center  of  Magnus  are  not  located  at  the  mass 
center. 
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The  terms  involving  accounts  for  the  center  of  pressure  being  located  oft  the  mass  center, 
white  the  terms  involving  accounts  for  the  center  of  Magnus  being  located  off  the  mass 
center,  (n  equations  I  and  2,  the  aerodynamic  coefficients  and  the  distances  from  the 
aerodynamic  force  components  to  the  projectile  mass  center  are  all  a  function  of  local  Mach 
number.  Typically,  in  flight  dynamic  trajectory  computer  codes,  this  dependence  on  Mach 
numbers  is  handled  through  a  table  look-up  scheme. 


4.  Aerodvnaidk  CoefTlclent  Lsiimation 
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The  lime- accurate  coupled  CFD/l^BD  simuiaiion  provides  a  full  tlow  solution,  including  the 
aerodynamic  portion  of  the  total  applied  force  and  moment  { X.Y,Z,L,M,N  ),  along  with  the 
tu  II  Slate  of  the  rigid  pro^ctlle  {  x,y,  z,  f.  6,  u,  v,  w,  p,  9,  r )  at  each  lime  step  in  the  solution 
The  ngid  state  of  the  prc^ectlle  Is  used  10  obtain  the  weight  portion  of  the  applied  force  so  that 
the  aerodynamic  force  can  be  isolated.  Using  the  information  provided  by  the  coupled 
CFDyRBD  simulation,  it  is  desired  to  compute  all  aerodynamic  coefTiclentsi  C,5,C,,.C^,C„<, 
^  3  n-stabi  lized  projecti  le,  the  Magnus  force  and  moment  are 

usually  sufficiently  small  so  that  C^,  and  C\«,  are  set  to  zero  and  removed  from  the  fining 
procedure  as  described  nexL 
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To  esDnuie  ihe  aerodynamic  coefTicienu  near  a  |>anicijlar  Mach  numbef,  a  set  of  rr  tune 
occuraie  coupleil  CFD/l^BD  sitnulaitons  are  creaieil  over  a  relatively  short  lime  period.  The 
initial  conditions  for  the  set  of  n  lime  histones  are  generated  lo  produce  a  rich  database  of 
aerodynamic  loads  and  projectile  stales  so  chat  a  unique  solution  can  be  obtained  for  the 
aerodynamic  coefficients.  The  initial  conditions  for  the  rigid  projectile  states  are  Ciaussian 
random  numbers,  with  a  mean  and  standard  deviation  selecied  to  cover  normal  operating 
conditions  for  the  projectile  Since  the  aerodynamic  coefficients  to  be  estimated  depend  on  local 
Mach  number,  the  set  of  ••  time  histones  is  repeated  at  m  differeni  Mach  numbers  of  interest. 
Thus  a  total  of  rr  short  time- accurate  coupled  CFD/RBD  uajectories  are  generated  to  support 

computation  of  a  complete  set  of  aerodynamic  coefficients  for  flight  dynamic  simulation. 

Since  initial  conditions  for  a  given  set  of  n  time  histones  are  randomly  generated  and  because 
Mach  number  changes  dunng  a  simulation,  Mach  number  varies  slightly  even  for  a  set  of  time 
histories  intended  to  be  generated  at  a  particular  Mach  number,  flence,  at  a  particular  Mach 
number,  all  aerodynamic  coefficients  are  assumed  to  vary  lineariy  with  Mach  number  For  a 
general  aerodynamic  coefficient,  this  variation  takes  the  form  shown  in  equation  Therefore 
parameters  at  intermediate  values  of  Mach  number  were  linearly  interpobted  as  shown, 

where  C  and  C,  are  the  aerodynamic  coefficient  values  at  Mach  numbers  slightly  (ess  than 
4  M  )  and  slightly  greatly  than  i  A/. )  the  target  Mach  number.  This  general  form  for  the 
aerodynamic  coefficients  is  then  substituted  into  the  aerodynamic  force  and  moment  equations 
Note  that  all  aerodynamic  coefficients  tliai  are  to  be  estimated  appear  in  the  force  and  momeni 
equations  in  a  linear  fashion,  suggesting  a  linear  least-squares  approach  to  estimate  ihe 
aerodynamic  coefficients  ai  each  Mach  number. 

Define  the  vectors  F^,  F^,  and  as  vectors  containing  all  the  unknown  aerodynamic 
coefficients  that  are  to  be  estimated  at  a  given  urges  Mach  number 
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Denote  the  total  number  of  unknowns  as  / .  For  a  fin-stabili2ed  projectile,  y  s  I4.  while  for  a 
spin-stabilaed  pro^ctile,  y  s  18 .  Assuming  each  lime  history  contains  k  time  simulation  output 
points,  then  A  "  n  linear  equations  in  j  unknowns  are  generated  at  each  target  Mach  number 
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Provided  ihe  matrices  A,,  ^ ,  A^,  and  are  maximal  rank,  a  unique  solution  for  ,  P,2> 

P^ ,  aivj  P^,  exists.  Thus,  propenies  of  these  matrices,  such  as  the  rank  or  singular  values,  can 
be  used  as  an  indicaior  of  ihe  suiiability  of  the  CFD/RBD  simulation  data  m  estimating  the 
aerodynamic  coefficients  at  the  target  Mach  number. 


S.  Kesulls 


In  order  to  exercise  the  previously  developed  method,  a  detailed  comparison  of  the  aerodynamic 
forces  and  moments  for  a  typical  finned  projectile  are  shown  in  figures  3-S  Mach  number 
varied  fiom  3.03  to  2.97,  covering  0.07  s  and  containing  3317  points.  The  projectile  had  the 
following  properties;  mass  of  4.  10  '  kg,  length  of  1.259  >  10  '  m,  diameter  of  I  319  ■  lO’ 

m,  and  axial  inertia  of  7.4  ■  10  ^  kg-m\  The  estimated  data  is  generated  using  equations  5  and  6 
with  the  identified  aerodynamic  coefficients.  Aerodynamic  forces  and  moments  agree  well, 
particularly  the  Y  and  Z  components,  which  appear  coincident  in  the  plotted  data.  The 
CFD/RBD  data  appears  slightly  noisy  in  both  the  axial  force  and  rolling  moment  However,  the 
estimated  data  removes  the  noise. 

Figures  9-17  present  estimation  results  for  an  example  of  a  finned  projectile  and  an  example  of  a 
spin-stabilized  projectile.  The  finned  projectile  is  a  120-mm  direci-fire  kinetic  energy  round, 
while  the  spin-stabilized  projectile  is  a  ISS-mm  shell  Synthetic  CFD/RBD  data  was  generated 
using  a  rigid  6  degree-of-freedom  trajectory  simulation.  The  solid  lutes  correspond  to  the 
coefficient  values  used  to  generate  the  synthetic  CFD/RBD  data.  The  “square”  symbols 
represent  the  finned  projectile  aerodynamic  coefficient  estimations,  and  the  “circle”  symbols 
represent  the  spin-stabilized  aerodynamic  coefficient  estimations.  In  figure  13,  Cmq  values  for 
the  finned  projeciile  are  scaled  by  a  factor  of  1/40,  so  that  data  for  both  rounds  is  easily  viewed. 
The  spin-siabilized  projectile  covers  a  Mach  range  from  0.6  to  4,  while  the  finned  projectile 
covers  a  Mach  range  from  I  to  4.  One-thousand  output  points  were  used  from  each  of  the  25 
runs  at  each  Mach  number  of  interest.  For  the  finned  projectile,  each  run  had  random  inidal 
conditions  for  angular  rates  with  a  zero  mean,  and  a  Ciaussian  standard  deviation  of  3  rad/s  for 
roll  rate  and  2  rad/s  for  pitch  and  yaw  rates.  For  the  spin-siabilized  projectile,  each  run  had 
random  iniiial  conditions  for  rol  I  rate,  with  a  mean  value  of  900  rad/s  and  Ciaussian  standard 
deviation  of  10  rad/s 
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Figure  3.  Estiftiaied  (dasfiirtf)  atMJ  CFU'RBD  (soUd)  body  sais  axial  force  vs.  time 


Figure  4  Eautitaied  (dashed)  and  CFD/RED  (seJid)  body  axis  side  tbree  vs.  utne 
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Figure  f>  Eeiimated  (dasbetf)  and  CFD.'HBD  (solid)  body  &iis  rolliog  momrm  vs  lime 


II 


Mz  (Nm) 


FigLK  7  Estinutcd  <das^c<f)  atMJ  CFB^HBD  ptuliiag  tnommi  vs  ome. 


FigLK  A  Estimated  (dasbetf)  artd  CFD.'HBD  (solid)  yawing  motneoi  ntne 
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Nole  *ipio*fUbiC)udpfDjeitile.  circle  •  peojeciJe.  $o3^  lioe  *  dei&«0\d  symboU  *catioule^ 

Figured.  OCOvs  Ma£l) nuftib^ 


t^oB  »quare  *^io«i£j^lBe^pr6^Dhe«  Ltfc3e  *  HrweJpeojeciile.  svlid  bos  *  ^u.  ef>d  symbols  ^BOmewd 
FigUM  10  CX2  VI  U^b  number 
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Mach  Number 


Nott  Ltfcle  *  Hrtf>cdpr4jKtile.  sdid  b 


M  *  wul  symbols  ^sDmoMd 


Figtif^  1 1  CNA  vs  ^tsdi  r^iimlwr 


Figiirs  13.  CLP  VI.  M^b  number. 
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Nott  »quars  *^io«i£j^lB^pr6^Dkf«  t9c)t  *  f^rwedprojeciile.  svlid  bos^  af>d  symbols  ^sDmsMd 

FigiMlJ  CLOD  Vi  Mi(^hDumbef 


Noto  square  *apio«d£jbilBO^pn^D^«  Ltf(3e  *  HrweJprojKble.  sdid  bos  *  af>d  symbols  ^somaMd 

FigiM  1 4  CUQ  VI.  Mdcb  number. 
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Ne4e  *  'prutiAbiliBd  pfcrMile,  arvK*  Afin^d  pre^eciilB. aelid l*ie  'tUa,  and  syinbc4>*  eiiifnAied 

Figure  IS  Suiion  liec  center  of  pmsuK  Mach  number 


Fifwe  16  CYPA  xi  Uecli  number  for  tpui-stabilized  pro|eciiic 
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Figure  IT  Sation  htK  cemer of  Uagttiks  force  vr  Mult  nuftibcr  for  sfHfi*siabi1i2Bd  projectile 


Paratnecric  trade  studies  were  conducted  to  determine  the  effect  of  the  number  of  runs  needed  for 
convergence  of  the  aeeodynamic  coefnciems,  os  well  as,  the  effect  of  the  length  of  the  time 
snippet  on  convergence.  Figures  1 1^23  present  results  for  CX2  and  CMQ  as  a  function  of  the 
number  of  runs  used  for  esiimodon.  The  number  of  runs  at  each  Mach  number  was  varied  to 
values  of  4,  6,  and  10  The  number  of  output  points  was  one  for  ell  runs.  Note  that  (or  10  runs, 
the  aerodynamic  coefficients  are  converged,  (n  figures  24  and  25,  the  effect  of  the  number  of 
data  points  in  each  time  snippet  is  investigated  Three  values  are  shown:  1, 2,  and  5  The 
number  of  runs  at  each  Mach  number  is  six,  and  the  snipet  length  es^uals  0  I  s.  Note  that  CX2 
and  CMQ  are  both  converged  with  five  data  points.  While  not  shown,  all  other  aerodynamic 
coefficients  are  converged  as  well. 
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Figure  19  CX2  V9.  Mach  numfin  ottedsia  poitn  pn  lime  anippet  aa  lime  snippets 
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FigLK  21  CMQ  U»ch  numiKr  onedaia  point  pn  onw  aniftpet  fourtinv  snippets 
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Cmq  I  Cmq 


FigiiK23.  CMQ  V9.  M act)  n umbel  ottedau  p«ijn  per  imc  sfuppet,  lObitie  snippet* 
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Figure  24.  CX2  Macb  Aufiibtf*  doia  pona  pet  linw  wuppct  ^  oot  rwo  (mddlcL  Hv^lnsM)^ 
irupp«u  s  1)^ 


FifiiM  2J  CUQ  vs  Macb  dumber  daia  poiiia  pcf  lime  saippet »  oee  rwe  (emddid)*  Ave  (ri^):  tma 
Aippeta  s  six. 


6.  Conclusions 


UsiBy  a  ume-acciiraie.  cotnpuuliona)  fluid  dyoamics  simubiion  tbai  is  iighity  coupled  lo  a  ngid 
body  dyniimics  sunuloiioo,  a  meibod  to  efDcimly  yeneraie  a  complaa  aerodynunic  description 
for  projectile  Qigbi  dyoamic  modeling  is  desenbed  A  set  of  n  very  ihon  une  snippets  of 
simulated  projectile  motion  at  m  difTereoi  Macb  nurebers  is  ci»npuied  aod  employed  as  baseline 
data  Tbe  comOiDcd  CFD/1U3D  analysis  computes  time  syfichroni?ed  air  loads  and  projectile 
Slate  vector  iiUbrreaiion,  leading  to  a  streigbeforuard  fining  procedure  to  obtain  tbe  aerodynasiie 
coefTicienis  Tbe  esisnation  procedure  decouples  into  four  subproblems  that  are  each  solved  via 
linear  least  squares.  By  inspecting  ibe  condition  number  of  each  fitting  matrix,  tbe  stolabtliiy  of 
tbe  lime  history  daia  to  predict  a  selected  set  of  aerodynamic  coefficients  can  be  assessed  Tbe 
overall  method  has  been  shown  to  work  well  for  both  fin-  aod  spm-siabiliaed  projeeiiks.  As 
would  be  expected,  convergence  of  the  aerodynamic  coefficients  is  strongly  influenced  by  the 
number  of  lime  snippets  and  ihe  number  of  points  in  each  time  snippet.  This  technique  provides 
a  new  means  for  the  CFD  analyst  to  predict  aerodynamic  coefTicdeols  for  flight  dynamic 
simulation  purposes.  While  CFD/RBD  is  computaiionaUy  intensive,  the  method  described  in 
this  report  provides  a  means  to  monitor  progress  in  csiimatmg  aerudynamic  coeflicienis  through 
tbe  individual  fittmg  matrices. 
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LUi  of  Symbols,  AbbreviarloDS,  and  Acronyms 


r.y.? 

U.VyW 


I' 

>' 

D 


a 

C,„ 

C,: 

C,, 


^•po 


•W 


Aalj 

CFD 

RANS 

RBD 

3-D 


Com[>onents  of  posioon  vector  of  tnajis  center  m  on  inertial  reference  frame 
Euler  roll,  pitch,  and  yaw  angles  of  box 

Components  of  velocity  vector  of  mass  center  m  body  reference  frame 
Components  of  angular  velocity  vector  in  body  reference  ftame 
Total  applied  force  components  in  body  reference  frame 

Total  applied  moment  components  about  mass  center  in  body  reference 
frame 

Magnitude  of  relative  aerodynamic  velocity  vector  of  mass  center 

Air  density 

Projectile  diameter 

Aerodynamic  angle  of  attack 

Zero  yaw  drag  aerodynamic  coelTicieni 

Yaw  drag  aerodynamic  coefficient 

Normal  force  due  to  angle  of  attack  aerodynamic  coefficient 

Magnus  force  aerodynamic  coefficient 

Roll  damping  aerodynamic  coefficient 

Fin  cant  aerodynamic  coefficient 

Pitch  damping  moment  aerodynamic  coefficient 

Distance  from  the  mass  center  to  the  center  of  pressure 

Distance  from  the  mass  center  to  the  center  of  Magnus 

Computational  fluid  dynamics 

Reynolds-averaged  Navier-Siokes 

Rigid  body  dynamic 

Three-dimensional 
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